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We use a functional approach to the Casimir effect in order to evaluate the exact vacuum energy 
for a real scalar field in d + 1 dimensions, in the presence of backgrounds that, in a particular 
limit, impose Dirichlet boundary conditions on one or two parallel surfaces. Outside of that limit, 
the backgrounds are described by a nonlocal effective action and may be thought of as modelling 
finite-width mirrors with frequency-dependent transmission and refiection coefficients. We obtain 
formal expressions for the Casimir energy in general backgrounds, and provide new explicit results 
in some particular cases. 

PACS numbers: 



I. INTRODUCTION 

The last years have seen a renewed interest in the Casimir effect tlf. The main reason for this has been 
the design and successful implementation of precision experiments The magnitudes involved in those 
experiments posed an important challenge to the first theoretical calculations, usually based on simplified 
theoretical models. Indeed, in those models, the materials were usually regarded as perfect mirrors, imposing 
strict boundary conditions on the quantum fields. 

In order to explain the experimental results, however, it is important to use more accurate models, including 
the corrections due to, for example, rugosity, and to finite temperature and conductivity. In this respect, 
there is a running controversy about how to model theoretically the finite conductivity of the material, and 
also about its influence on the Casimir forces Q. The computation of the self-energies due to quantum 
fluctuations is even more controversial: while the Casimir forces between two objects become finite and 
independent on the properties of the materials in the limit of perfect conductivity, the situation for the 
self-energy is, by far, not so clear. The reason may be traced back to the fact that the local energy density 
diverges close to an idealized boundary and it is therefore necessary to introduce surface counterterms 
in order to get finite results for the self-energy Since the self-energy is relevant, for instance, to analyze 
physical problems like gravitational effects produced by the vacuum energy, a careful analysis of the sources 
of divergences in the stress tensor as well as in the total energy is in order. 

Motivated by those problems, in this paper we are concerned with the calculation of the vacuum energy 
distortion that results from the presence of some particular backgrounds, which are meant to represent more 
realistic boundaries, for a real scalar field theory in d + 1 spacetime dimensions. Following @, we will 
choose backgrounds which may be used to reproduce exact Dirichlet boundary conditions on flat surfaces, 
when a particular limit for the parameters that define the background is taken (see also 0]). As in 
we also introduce two properties in the descriptions of the plates: finite widths and finite strengths. The 
former is represented by a smooth function a^, depending on a single variable (the normal coordinate to 
the plate), which has a characteristic width 2e. Regarding the strength of the coupling between the walls 
and the fluctuating field, we extend the kind of background considered in by allowing for a frequency 
and (parallel) momentum dependent form for a function A, which determines the strength of that coupling. 
Any non-trivial dependence on those variables shall amount, in coordinate space, to a nonlocal form for the 
action that describes the interaction between field and plates. The main new contribution of this paper is the 
derivation of expressions for the Casimir energy, taking into account nonlocal effective actions that model 
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finite-width mirrors with frequency dependent reflection and transmission coefficients, using the functional 
formalism as a tool. We do that for general cases, obtaining formal expressions, as well as for particular 
examples, where we derive more concrete results. 

We will also be able to consider smooth backgrounds, in order to test the dependence of the Casimir 
energy with the 'sharpness' of the boundaries. As already stressed, these generalizations are of interest for 
the calculation of the Casimir forces in realistic situations, and also from a formal point of view, since many of 
the divergences that appear in the calculations of the zero point self-energies are related to the introduction of 
sharp interfaces and/or ideal boundary conditions. In this way, with finite-widths and frequency- momentum 
dependent strengths, one can simultaneously control the two sources of UV divergences that pop-up during 
the calculation of vacuum energies. This will also allow us to avoid the unnecessary introduction of ill-defined 
quantities at the intermediate steps, in the course of calculating physical observables. 

The calculation of the Casimir energy in the presence of realistic mirrors is a subject that have been con- 
sidered in a large number of previous papers. In particular, within the Lifshitz approach it is. in principle, 
possible to compute the zero point energies for slabs of arbitrary width and arbitrary electromagnetic prop- 
erties, bounded by flat interfaces 0]. In our approach, we will be able to consider smooth surfaces, including 
some which approximate slabs (to an arbitrary degree). Moreover, the functional method is particularly 
well adapted to consider the effects on the zero point energy of the degrees of freedom inside the material. 
Indeed, the nonlocal effective action may be thought of as resulting precisely from the integration of those 
degrees of freedom. Thus our initial nonlocal effective action can be considered as a toy model for the 
interaction of the would-be electromagnetic field (here the quantum scalar field) and the charges confined 
in the material. In order to illustrate the method, in the present paper we will consider a particular class 
of nonlocal effective actions for a scalar field, leaving the generalization to the electromagnetic case and the 
derivation of the effective action from first principles for a future work. Nonlocal effective actions have also 
been briefiy considered by other authors, see for instance 0]. 

^From a more technical point of view, we shall use an extended version of a previously used functional 
approach Q to the Casimir effect with perfect mirrors, in order to cope with the more general situations 
considered here. 

The organization of this paper is as follows: in SectionlTTl we introduce the general method used to calculate 
Casimir energies. We obtain formal expressions for the total vacuum energy for one and two mirrors. We 
then present explicit results: in Section IlIIl for the change in the vacuum energy under the introduction of 
a single mirror, to consider afterwards, in Section llVl the case of two mirrors. For the single mirror case, 
we compute the Casimir energy associated to smooth and piecewise constant backgrounds. We discuss the 
dependence of the energy with the width e, the cutoff frequency, and the sharpness of the boundary. For the 
case of two mirrors, we present a general expression, suitable for numerical evaluation, and derive some of 
its general properties. The expression becomes much more simple for the particular situation in which the 
width of the mirrors is much smaller than the distance between them, but not necessarily smaller than the 
inverse of the strength of the coupling A^^. In Section |V] we summarize our conclusions. 

II. THE METHOD 

In this Section we set up the problem, and derive the general expressions subsequently used for the 
calculation of vacuum energies in some specific cases. 

We shall consider either one or two identical, fiat, parallel and finite width (~ 2e) mirrors in d spatial 
dimensions. A coordinate system has been chosen such that their 'centers' correspond to = for the case 
of a single mirror, and to = and Xd = a, for the case of two parallel mirrors. In both cases, the system 
will be described by an Euclidean action: 



(1) 
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where 5*0 defines the free theory: 



Soiv) = I I d^+'x d^^d^^ , (2) 



while Si{ip) is a term, quadratic in the field, that contains the interaction with the walls. Introducing a 
number A'^ = 1, 2 depending on whether one has one or two walls, respectively, we have: 



N 



^/ = E^i"^ (3) 



where: 



S\°'\f) = '^ J dxo J dx'„ j d'^ ^x\\ j d'^ '^x\^ j dx^ ^'(xo, X|| , x^) A(xo-a;o; X|| -xJ|)CT£(xd-aQ) v3(a;o, xj| , x^) , 

(4) 

with fli = and 02 = a. Xj| denotes the 1 coordinates parallel to the mirror: Xi, 0:2, ■ ■ • , Xd-i- We regard 
the interaction term as an effective action coming from the integration of the degrees of freedom confined to 
the walls that interact with the scalar field 1^3(2;). On general grounds, the quadratic part of that effective 
action should be of the form: 

5cff = i y d'^+^x j d'^+'x' ^ix)T2{x;x'Mx') . (5) 

Taking into account translation invariance in xq and X|| , the dependence of the kernel T2 on the coordinates 
must be of the form: 

T2ix;x') = T2ixd,x'd,xo - Xo,X|| - x'||) . (6) 

If, in addition, we assume locality in the perpendicular coordinate Xd, we arrive at ([4]), where the functions 
A and Ce encode the structure of the walls, such as their reflection and transmission coefficients and their 
widths. We will assume, without any lose of generality, that A is an even function (any odd part would cancel 
away in the action). Therefore its Fourier transform. A, is real. Besides, we will also assume that A is strictly 
positive, so that the interaction with the walls is repulsive at all frequencies. Note that the nonlocality in 
configuration space will produce reflection and transmission coefficients dependent on the frequency uj, as 
well as and on the wave vector along the parallel plane, ky . 

On the other hand, cr^ is an even, strictly positive and continuous function, approximately constant on a 
region of size ^ 2e around 0. There is no big loose of generality by assuming strict positivity for this function. 
Indeed, assuming that one wanted to consider, for example, a function whose support is the interval [— e,e], 
we could approximate it by a strictly positive one that vanished, arbitrarily fast, outside that interval (see 
Section IlIII where this case is dealt with in some detail). Since appears multiplied by A, we may impose 
the condition: 

+00 

dxd(7c{xd) = 1 , (7) 

-00 

without restricting the actual interaction at all. Rather, it is a way of disentangling the 'shape' (attributed 
to (Tc) from the strength (carried by A) of the interaction. 

In our conventions, the Euclidean coordinates shall be denoted as x^, /x = 0, 1, . . . , d. The Xd coordinate 
points along the normal direction to the walls, while X|| = (xi, . . . , Xd-i) are parallel to them. The d spatial 
coordinates are collectively denoted by x. 

We now consider the vacuum energy cost, of distorting the vacuum by the introduction of the mirrors. 
This quantity may be written as follows: 
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where 

Z = Jvipe-^^'^^ , Zo= Jvipe-^"^'P\ (9) 

and T is the extension of the (imaginary) time interval. 

In more than one spatial dimension, Eq is proportional to the 'area ' L'^-i of the mirrors (assumed to be 
{d — l)-dimcnsional squares of side L). Since L ^ oo, it is convenient to introduce the energy density Eq, 
such that 



£o = lim ln(— ). (10) 



Besides, in the case of two mirrors, one is usually interested not in but rather in a subtracted quantity, 
£0, defined as the difference: 

£0 = £0 - £o{^) (11) 

where £q{oo) denotes the surface energy density when the mirrors are separated by an infinite distance. 

£q is finite even if ideal mirrors (i.e., imposing Dirichlct boundary conditions) were considered. The 
reason is that the self-energies of the mirrors, being translation invariant, are cancelled when subtracting 
the a ^ 00 contribution. This, however, would involve a regulatization in order to avoid the having to deal 
with the difference between two ill-defined (divergent) quantities. In our case, the self-energies are finite due 
to the presence of a physical rcgularization mechanism. Therefore, the subtraction of the self-energies is a 
well-defined step, without having to invoke any additional regulator. 

To proceed, we note that a quite natural extension of the functional approach followed in previous works 
can be implemented here. Indeed, we can introduce an auxiliary field, ^a, in order to rewrite the exponential 
factor S^"\^): 

g-S<°'(>^) ^ _1 J Jy^^^~^fd''+'xfd''+'x'i^ix)X-\xo-x'„■,^^^-^\^)cr.{xa-a„)S{x,-x'^)i^(x') 

X I d''^^^io,{x)a-^{xd-a^)ip{x) ^ ^-[^2) 

where the factor A/", independent of a, is given by: 

J\f ^ J j)^^^-^ f d'^'^^x f d''+^x'i^{x)X-^{xo-x'g-xii-x\^)a,{xa)S{xa-x'j^)ic,{x') ^ ^j^g-j 

and we have used A~^(a;o — Xq; X|| — x'|) as a notation for the inverse kernel associated to A(xo — x'^; X|| — xj| ), 



I.e. 



A (.To 



{,';x||-x|pA-i(x;;-4;xjf-xI|) = / A-i(a;o - xf,'; x,, - xj[)A(a;;,' - xf,; xj[ - xj,) 



.11 II 

'-^11 " ■^o '-^11 



= (5(a-o-x[,)J(x|| -xj|) . (14) 

An important difference with the case of ideal, zero-width mirrors, is in that the auxiliary fields now live 
in d + 1 dimensions, rather than on the d-dimensional submanifold = 0. The present case reduces to the 
ideal situation if is replaced by a (5 function (of which it may be regarded as an approximation when e 
is finite). In other words, there is a 'dimensional reduction' in the auxiliary fields when e — > 0. Particular 
profiles for the function are introduced in Sections IIIII and IIVI obtaining for them explicit results. 

We now use ([T^ to rewrite Z and afterwards integrate out the scalar field (p, to obtain: 

1 r ^ 

Z ^ Z^ j Y[V^^ e-Hrf^+^-/rf^+^-'E^,,., , (15) 
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with a matrix kernel whose elements flaf3 are defined by: 

+ cr^{xd - Ua) A{xo, :s.i\, Xd; x'q, x' \i, x'd)aeix'd - afs) , (16) 
where A is the free scalar-field propagator: 

We then make, for each auxiliary field ^q,, the redefinition: 



^a{xo,xi^,Xd) ^ / A2 (xo - a;o;x|| - xJ|)^a(a;o,x'|,Xd - aa) , (18) 

both in the explicit integral over ^ in (fT5|) and in the implicit one in the definition of Af. We have assumed 
that the square-root kernel is well-defined, what is consistent with the assumptions about its properties. 
This redefinition produces a non-trivial (A-dependent) Jacobian for each factor P^q,. An identical Jacobian 
does, however, appear also in each denominator Af, and therefore they cancel each other. Besides, after the 
redefinition, the J\f factor changes: Af Af, which is independent of A. Thus, after a trivial shift in the 
coordinates, 

{Af)^ J 

where: 

n^pix,x') = da.[3 S''''+^\x - x')a,{xd) 

+ <7c{xd) / A^(2:o -a;Q;x|| -xJ|)A(a;o,x'[, Xo",xJ|',a:;d + a^j) 

X Xi{x'^'^x',;x\l'-x\^)a4x'd). (20) 

We then proceed to perform yet another redefinition of the auxiliary fields, now involving a diffeomorphism 
in the Xd coordinate. We introduce the one-to-one mapping Xd ^ z that results from the differential equation: 

= 'Jeixd) , (21) 

dxd 

where the assumed positivity of ctc comes in handy to ensure that the change of variables is non-singular. 
Imposing the condition z(0) = 0, the solution to the previous equation is imique. We also note that, as a 
consequence of the change of variables, 

S{xd-x'd) = cr^{xd) S{z - z') , dxd(T^{xd) = dz . (22) 

Equation (|21[) . together with the condition z(0) = and the normalization condition for CTj imply that the 
range of z is (regardless of the particular profile used for a^) the finite interval [—5, 

Keeping the same notation for the auxiliary fields when written in terms of the new variables, wc then 
have, in a simplified notation: 

Z = Zo [ f[l^^ g-3/xo.-||,-;xi,xJ|,.'5:<,,3«o(:ro,X|,,z)K:„3(a:^o,X||,^;a:^;„xj|,^')«,3K,x;|,2') ^^^^ 
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with 

/Ca;3(a;o,X||, z;xo,x'|, z') = (5^/3 (5(xo - a;o)5(''~^)(x|| - x'|)(5(z - z') 

X^xq - a;o;x|| - x'[)A(xo , X|| , a;rf(z) + tta] Xq", x'| , a;^(z') + a^)A5(a::o" - a;o;xJ|' - X||](24) 



We note that, had any Jacobian arisen because of the last field diffeomorphism, it would have, again, been 
cancelled against an equal object coming from JV. Furthermore, since after the redefinition Af became 
independent of A and cr^, that factor has been dropped. 

Then we Fourier transform all the spacetime coordinates for which there is translation invariance, namely, 
a;o,X||. The transformed operator has the following form: 

i^ai3{z;z') = SapS{z - z') + X{uj,ki\) Dai3{z;z') , (25) 

where 

A(w,k||)= dxa I d'^-i.Tii e-<'""''+''ii"ii)A(xo,X||) (26) 

J —oo J 

and 

-K\xd(z)-Xd(z')+aa-afi\ 

DMz; z') ^ , (27) 

where k = ^J^"^ + ^j^ (to simplify the notation, we omit writing the dependence of D^piz; z') on fc|| and uj 
explicitly) . 

Then we have the following expression for Z: 

Z = Zq (det^)"^ , (28) 

and: 

So = lim 3-^Tr(lnK;), (29) 

where 'Tr' denotes trace over all the variables (including uj and fc||). Introducing the symbol 'Tr' for the 
(reduced) trace over the Hilbert space of functions depending on 2 G [— i, i], wc may write a more explicit 
formula, where the trace over the variables for which the kernel is translation-invariant is explicit: 

^o^y^y(^Tr(ln/C). (30) 

/,From the last expression one can, in principle, extract the vacuum energy relevant to each case. However, 
since there are important differences between them, we present now separate calculations of corresponding 
to two different physical situations: for the case of one mirror, and for a two-mirror system. 

A. One mirror 

This case amounts to a single index a = 1; thus we suppress the a, /3 indices, and consider just the kernel: 

K.{z;z') = 8{z- z') ^\(uj,\\)D{z,z') , (31) 
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with 

D{z;z')^ . (32) 

We shall calculate the trace of the logarithm of JC, by finding the eigenvalues of /C. The equation for tpa{z), 
the eigenfunction associated to the eigenvalue a, may be written as follows: 

' dz'X{uj,kn)D{z;z')'iPa{z') = (a - IjV'al^) • (33) 

This can be converted into a differential equation, taking into account the fact that D{z; z') is obtained by 
performing a Fourier transformation plus a change of variables in the Green's function for the real scalar 
field. Indeed, starting from: 

( ~ ^ + '*^) ^o(a;d, x'^) = S{xd - x'^) , (34) 
where Goixd,x'j) is the Fourier transformed of the free propagator: 

Go{xd,x'j^) = — , (35) 

we obtain, by changing variables: 

[-ae(z)|-(a,(z)|-) + ^]Diz:z') = d^z) S{z - z') , (36) 



dz dz 

where 



ae{z) = ae{xd{z)) . (37) 



We see that, by acting with the differential operator that appears on the left hand side of ([36|) on both 
sides of p3p . it becomes a differential equation with the structure: 



L[i;^]{z) = , (38) 

where L denotes a linear differential operator of the Sturm-Liouville type: 

L[f]{z) = -^[pW^] + : (39) 

with 

^2 

p{z) = a,{z) , q{z) = ^ , (40) 

and the eigenvalue of L, ^a, determined a through the relation ^q, = ^jzrj-- The eigensystem corresponding to 
the operator whose determinant we need has the coefficient functions p and q, which are in turn determined 
by the properties of the mirror. 

An important fact to note is that the space of functions to be used for the calculation of the eigenvalues 
of L corresponds to functions that vanish at 2; = ±i. Indeed, the K operator is self-adjoint under the L^(R) 
scalar product defined by: 

r+co 

if,g) = \ dxdr{xd)a(xd) (41) 
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for any pair of square integrable functions: /, g. This scalar product becomes, when written in terms of the 
new variable z: 

^ / ^^nz)g{z). (42) 

In particular, this implies that normalizablc functions in the interval [— ^, \] must vanish when z = ±1/2 
(of course, they should also decrease sufficiently fast for the integral to converge). The reason is that, from 
the defining properties of ct^, it follows that vanishes at least linearly at those points. 
It also implies that eigenfunctions of the Sturm-Liouville operator corresponding to different eigenvalues 
shall be orthogonal, for the scalar product ([42]) . which includes the weight function I/cte. 

Thus, we conclude that, to evaluate the vacuum energy in the one mirror case, this procedure has lead us 
to the expression 

''■ = \ltl$^T.H'n(^-h.4^ (43) 

where / labels the eigenvalues a. They are, in general, non trivial functions of their arguments, and they 
depend on the particular function considered. 



B. Two mirrors 



We shall now consider the situation of two identical mirrors, located at a distance a apart. Thus, we are 
concerned with the quantity: 

£o(a) = fo(a) - fo(«3) (44) 
We have to consider now the 2x2 matrix kernel: 

where the two diagonal elements in the expression above coincide with the (identically noted) kernel corre- 
sponding to the single- wall case, equation ([3T|l . There appear also two new kernels V and Q: 



^-K.\xd(z)-Xd(z')-a\ 

V{z-z') ^ A(cj,k||)go(:crf(z) + 0;xd(z') + a) = A(w,k||) — , (46) 

and Q = V\a^-a- 

One can then use the expression for the determinant of a matrix in terms of a combination of its blocks, 

Aet[K,j\ = (det^)^ det(l-O) (47) 

where 

O = iC-^ViC-^Q . (48) 

Then we note that, since det/C is independent of a, the Casimir energy density, obtained by subtracting the 
contribution with the plates at an infinite distance, becomes: 

£o(a)= lim^^Trln(l-O). (49) 
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It is important to note that, in the last equation, the a-independent contribution, corresponding to the 
self-energies of the plates, have been completely subtracted. Note that that contribution yields exactly twice 
the vacuum energy for one plate. Of course, had A been assumed to be a constant, those self energies would 
have been, generally, divergent. 

Thus the UV behaviour of is milder that in other approaches, where those contributions have to be 
dealt with in order to give sense to the Casimir energy. Here, as wc shall see in the examples, it is already 
finite. Assuming that fc|| , e, a) are the eigenvalues of O: 

'^^'^ -ir^i ^Ei-[i-^'(-'^ih^'-)] ' (50) 

the expression that we will use as a starting point for the calculation of the Casimir energy density in concrete 
examples. 

The finiteness of ([50]) will follow from the fact that, as a function of the frequency and momenta, O is 
bounded by the exponentially decreasing factors of those variables, carried by V and Q. 



III. RESULTS FOR THE TOTAL VACUUM ENERGY IN THE SINGLE-MIRROR CASE 



A. Smooth (Te 



Wc first assume a particular form for the function a^- A very convenient choice, from the calculational 
point of view is: 

f^eixd) = — sech^( — ), (51) 
2e e 

where the ^ factor has been introduced in order to satisfy dxd(J^{xd) = 1, the normalization condition 
corresponding to an approximant of the (5-function. Of course, the wall is then localized around Xd = 0, with 
a width It. Besides, for this profile we may find the explicit relation between Xd and z: 

Xd{z) — e arctanh(2z) 

z = itanh(^) . (52) 



Note that z S [— ^, We immediately find: 

1 



a,(z) = ^[l-(2^)'] • (53) 



Applying the differential operator on the Ihs of ((36|) to both sides of ((33)) . and then changing variables: 
z ^ u= 2^, we obtain a differential equation for the eigenfunctions: 

-[(l-.^)-^.(.)] + [^(Hy-H^]^.H - 0. (54) 

We recognize the associated Legendre equation, which has regular independent solutions when: 

e^(^,k||) 

^^^^(^ + 1), ^ = 0,1,... 

CK = m, m=l,...,Z. (55) 
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Note that the m = solutions have been discarded, since the last condition does not lead to any non-trivial 
solution. The eigenvectors are then the polynomials P™(2z) which (since m ^ 0) vanish at z = ±i, as 
expected from our general analysis at the end of the previous section. 
^From the above, we conclude that the eigenvalues are: 

a = + (56) 



while, for each I, we have the constraints: 



fc2 = l,...,Z, (57) 



which restrict the allowed values of fc|| = |k||| and lo. We may then use the explicit form of the eigenvalues 
to obtain: 



Jd-l I 



(=1 m=l J \ ) 

X ^e^^^-H In [1 + ^1^], (58) 

Note that we are writing a general expression as a function of the dimension of space, d, but we are not 
meaning by that that a dimensional regularization is used. Indeed, as we will see, this expression will be 
regularized by the function A. 

In the case d = 1, there is no integral over fc||. Furthermore, the integral over w becomes trivial because 
of the constraint, and thus we derive an expression for the vacuum energy in terms of a double sum: 



-.00/ 



;=i m=l 



^A(f) 
2/(/ + 1) 

6A(-f 



or, 



1=1 m=l ^ ' 

since A(a;o — Xq) is assumed to be real (i.e., no dissipative coupling to the mirror). 

The result lpD|) shows an interesting relationship between the large-w behaviour of A and the size of the 
mirror, e. Indeed, assume that A becomes negligible small above some cutoff frequency Wc, then we have the 
approximate expression: 

00 min{/,[ewe]} \(m\ 
(=1 ?ri=l ^ ' 

Note that the energy is zero if e < l/wc. This may be understood from the fact that no modes could be 
trapped inside the mirror, since it is transparent for wavelengths smaller than e. 

To proceed in the case d > 1, we will assume that A(ti',k||) = A(k). The motivation for this assumption is 
that, with this choice, the reflection and transmission coefficients will depend on the wave number fc^ in the 
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perpendicular direction (note that when rotated back to Minkowski space, k becomes kd)- In this case, the 
integrals in Ea. ((58)l can be easily computed. We obtain 

^ ' 1=1 m=l ^ ' 

where Sl^-i = ••• is the solid angle in d — 1 dimensions. If there is a cutoff Wc, as in the previous example, 
and besides we have that for k < o^c it is constant: A = Aq, then: 

Of course, we also have the same relationship between Wc and e as in the previous example with = 1, since 
the defect is essentially one-dimcnsional. 

Let us see that the introduction of the cutoff lo^ produces a finite result, in any number of dimensions. To 
see this, we just consider the case of a A which is constant below the cutoff, and zero above. Then we may 
rewrite (|63p more explicitly, as follows: 

+ y inri + _!^l. (64) 

2(27re)'*^^ ^ ^ 21(1 + \y ^ ' 

Now, the first term on the rhs above is a finite sum, while the second one involves a convergent series, so 
the result is finite. Of course, the same holds true for the case of a (continuous) A function which is not 
necessarily constant below the cutoff. Indeed, one may simple bound that function by a constant, and then 
use the previous result. 

B. Piecewise constant profile 

We consider now another important profile for the function a^. Wc assume that this function takes the 
approximately constant value l/(2e) inside the interval [— ^, 5], and vanishes very fast outside. The function 
is then extremely smooth, except for small intervals around x^. = ±e where all the variation is concentrated. 
Of course, one can approximate a 'square barrier' function with this kind of profile, and this is the motivation 
for considering it. 

In the limit when that profile is approached, we are left with the equation for the eigenvalues; 

- ^M^) + (26«)Va(z) = ^^^^^o (65) 

in the interval (— ^, ^). With the boundary conditions we find the eigenvalues of K., which form a discrete 
set: 

2eA(cj,k||) , , 



Thus, in d spatial dimensions, the expression for £ becomes: 



^ _ 1 [ duj [ d''-^k\\ ^ 
^° - 2 i 2^7 (2^ A. 



2eA(w,k||) 



P7r2 + (2e)2(tj2 + fc2) 



(67) 
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which differs from the one obtained for the other profile used for ctj. It sliould be noted that the procedure 
above for the square barrier profile relies upon the existence of a region where that function is approximately 
constant. This would not make sense if we wanted to consider the i5-function (e — > 0) limit, which we shall 
study (sec next section) using the (sech)^ function instead. 

We assume, as stated before, that A depends on lo and fc|| trough the particular combination k = 
i.e., A — X{k). Then we may write a more explicit formula for £q: 



II' 



£o 



2d^d/2r(d/2) 



dKK'^-^ In 
1=1 



2e~X{K) 



(68) 



The series is convergent, and one can immediately find its large-K behaviour: 



1=1 



In 



2eA(K) 



P7r2 + (2eK)2 



A(^) 



(69) 



^From the previous expression, we see that if A vanishes above a cutoff, the energy is finite. We see that 
it would also be finite if it vanished following an exponential law, or even a power-law: A ~ with 
a > d - 1. 



C. UV divergences 

We conclude this section by presenting a brief study of the UV behaviour (divergences) of the expressions 
that we used to calculate the energies in the one mirror case; i.e., of the self-energies, when one approaches 
the constant-A case. 

Wc will study here that limit starting form the self-energy for a finite-width mirror with a cutoff for k, so 
that the energy is finite, and then take the relevant limit, which corresponds to a constant A (that eventually 
tends to infinity), and a vanishing e. The most transparent way to take these two limits is by letting first 
e — > 0, and then making A go to a constant. The advantage of using this order is that the first limit can be 
taken exactly. Indeed, when e — > 0, we have: 

IC{z;z') ^ 6{z^z') + ^ , (70) 



or, in operatorial form. 



^ - (1 + + V , (71) 



where p and i] are projection operators whose kernels are: 

piz;z') = 1 , 77(z;z') = S{z - z') ~ p^z, z') . (72) 
Thus, the expression for the vacuum energy in this limit is 

\£o] , = , , , , r dnn"-^ Inh+Ml . (73) 
L uj.^o 2'^7r''/2r(d/2) io ^ 2k ^ ^ ' 

Note that the large-K behaviour of the integral is: 

' dnn"-^ In [1 + 4^] ^ \ H dn k^-^~X{i^) , (74) 

2k 2 ./n 
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as in the sharp-cr case. Assuming a constant A, and introducing an Euclidean cutoff A for the K-integration, 
we see that: 

' dn k"-' In [1 + ^1 ^ -Xo r dn n"-^ , (75) 

where the ^ refers to the large-K behaviour of the integrals only. The UV divergences in the previous 
expression are then quite immediate to extract: in = 1 (two spacetime dimensions) there is a logarithmic 
divergence: 



where is a momentum scale. It should be noted that these divergences may be interpreted in quantum field 
theoretic terms. Indeed, when A is a constant, the model we consider may be regarded as a free scalar field 
theory with Lp^ insertions; the latter coming from an expansion of the term proportional to A in the action 
(which is local when A is a constant). For n insertions (i.e., the term of order A") in d + 1 dimensions, the 
superficial degree of divergence (5 of a term contributing to 8q with n insertions of (f^ is, by power counting: 

5^{d+l) + n{d-l)-n{d+l) . (77) 

For the d = 1 case, we see that only the 7i = 1 term (linear in A) diverges and it does so logarithmically, as 
we have seen above. In two spatial dimensions, there is only a linear divergence, in the n = 1 term, while 
for d = 3 there are divergences for n — I (quadratic) and n — 2 (logarithmic). Therefore, we expect to have 
one, two and three renormalization conditions for d = 1, d = 2 and d = 3, respectively. 

To make sense of those divergences, rather that dealing with the one particle irreducible functions corre- 
sponding to the operator insertions, we apply the renormalization program in a more straightforward way, 
in terms of the vacuum energy, as follows: subtracting from the integrand in the expression for £q with a 
constant A its MacLaurin expansion in A up to an order 5+1, where S corresponds to the higher degree of 
divergence at the given d, we have a convergent expression. 

The renormalized energy will then have a finite degree polynomial, whose coefficients must be fixed by 
imposing suitable renormalization conditions. For example, in c? = 1: 

£o = ci Ao + A£o , (78) 

where ci is a constant (depending logarithmically on the cutoff), and ASq is the finite expression: 

Ao \ Ao " 



1 

A£o ^ -z- I dn 
271- Jo 



^ 2k/ 2k 



(79) 



A simple rescaling shows that A£q is also linear in Ao (although with a cutoff- independent, finite coefficient). 
Thus we conclude that the renormalized energy will have the form: 

[^"]d=2 = ^1^0 ' (80) 
where Ci requires, to be fixed, to know the vacuum energy at some scale. In ci = 2, on the other hand, 

£o = ci Ao + C2 A^ + A£o , (81) 

where now, 



1 r°° 

A£o = — / dKK 
47r Jo 
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which goes hke Aq; then: 

[£o]_ = CiAo + C2A^ (83) 
Finally, in d = 3, a similar procedure yields: 

[fo]_ = C1A0 + C2A2+C3A3. (84) 

As already mentioned, the vacuum energy is simply related to the effective action in the presence of (ys^ 
operator insertions; hence the renormalization conditions on the vacuum energy may be related to conditions 
for the corresponding one particle irreducible functions. From a more phenomenological point of view, the 
meaning of the renormalization conditions is best put in a negative way: they summarize the ignorance one 
has about the vacuum energy, when it is divergent. Namely, one can predict the dependence of the vacuum 
energy with the coupling constant, except for the first few terms in a McLaurin expansion of the energy 
in terms of Ao, the k 00 part of that coupling constant. It is perhaps worth emphasizing that in a real 
situation there must be a cutoff; hence Aq = 0, and the divergences above are replaced by finite, cutoff 
dependent terms, which one does not need to renormalize, and whose precise form depends on the details 
on the defect (its profile, for example). 

We sec that wc would need to impose three of those conditions in order to completely fix the renormalization 
constants C,;. We agree with the results of 0], for the case of a sharp defect. 

Finally, note that the would be n = (A-independent) divergences do not appear because we measure 
energies with respect to the vacuum in the absence of mirrors. 



IV. RESULTS ABOUT THE CASIMIR ENERGY FOR TWO MIRRORS 



We shall now consider the evaluation of the Casimir energy, as a function of the different parameters, for 
different profiles. 



A. Thin mirrors 



As a first check of the method, we consider firstly the case of thin walls (e 0) and arbitrary A. This 
corresponds, physically, to e much smaller than the other two parameters with the dimensions of a length, 
a and A~^. Later on we shall also impose the condition that A tends to infinity to recover the well-known 
Dirichlet case. 

The e — > condition can be easily imposed on /C; indeed, when e — > 0, one has the following expression 
for that kernel: 

JC-^l,,/)^Siz-z')--^H— (85) 

2k + \[K) 



and 



nz;z') - ^e--^. (86) 



2k 



Then it is immediate to see that, in the same limit, 



0{z;z') - ( (87) 
2k + X[K) ' 



15 



i.e., it becomes independent of z and z' . Thus the Casimir energy density for e — *■ is given by the expression; 



£o(a) 



1 



2^7r^/2r(d/2) 7o 



dKK'^-^ In 



1 - 



'2h. + A(k)' 



\2 g-2rea 



which could be evaluated numerically for different interesting functional forms of A, depending on the material 
considered. 

We can obtain exact results using some particular profiles. The Dirichlet case is obtained by considering a 
constant A which tends to infinity. Note that the integral is convergent for any finite A, thus the limit could 
be taken after evaluating the integral over momenta. Nevertheless, taking the A — > c» limit before we can 
perform the integral exactly, in any number of dimensions, d. For example. 



^o(a) = 



C 



167ra2 
" 1440ai' 



for = 1 
for = 2 
for d = 3 



(89) 



There is another profile for A which allows we to find exact results, albeit it is unphysical regarding the 
properties of the function A. However, it serves the purpose of illustrating a property of the Casimir effect in 
the Dirichlet case: consider the profile A(k) = la k, where a is a constant, which has the unphysical property 
of growing with the frequency and momentum. In this case we have: 



1 



I{d,a) , 



dK K 



d-1 



In 



1 - 



(90) 



where 



I{d,a) 



2dn'i/2Y{d/2) Jo 



dxx'^-^ In 



(91) 



is a finite number, depending on the constant a and the dimension, d. Note that this kind of profile produces 
a dependence of the energy with the distance that is identical to the Dirichlet case, although with a smaller 
coefiicient, in spite of the fact that the coupling constant is a function that grows with k. This is a reflection of 
the fact that, although A is not infinite, its particular form introduces reflection and transmission coefficients 
that are independent of k, as in the Dirichlet case. 

Of course, when a oo, we recover the Dirichlet result, which shows that one can approach it not just 
from the constant A case, but also starting from a rather different profile. 

If, on the other hand, e is negligible with respect to a, but not necessarily in comparison with A^^, we 
may use the approximation: 

\xd{z) ~ Xd[z') + a\ - a (92) 
in the expressions defining V and Q. Then, from the definition of O, we see that: 



0(^;2')-^(^) [ dz2 dz3K:-\z2;z3)\ dz^JC-Hz; zi) e 



which is independent of z' . 

Then, O has only one non-vanishing eigenvalue, i?, for each k, as in the e 



case: 



■d 



V 2h 



dz / dz'K. ^{z;z') 



-2Ka 



(94) 
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and 

Of course, the result depends on e because of the object J dz J dz'IC^^{z; z'). For the case of the piecewise 
constant defect, for example, we find: 

fdzfdz'IC~\z;z') . ^V,/ , (2fc + l)^.^ + (2e.)^^ 
J J TT ^ (2fc + 1)2 (2fc + 1)2^2 + (2eK)2 + 2eA(«;) 

^ ^ (2/c + 1)2 (2A: + 1)2^2 + (2eK)2 + 2eA(K) ' 

We conclude the study of the thin wall case by mentioning that it is possible to study analytically the 
nct-to-leading term in an expansion in powers of e for £0(^)1 by extending the method used in thin- wall 
approximation. A rather lengthy calculation shows that the first-order term exactly vanishes. Thus, the first 
non-trivial correction to £o{a) (if present), can be of the form Ce^/a'^'^'^, in d spatial dimensions. 

We also wish to pinpoint to a phenomenon, which appears when considering the finite size limit in expres- 
sions obtained by using the eigenvalues corresponding to a particular profile. The e — > and A ^ 00 limits, 
which we have studied exactly at the beginning of llV Al have been found by using the explicit form of /C~^. 
It corresponds, in fact, to the exact solution when the mirrors are J-like potentials. That limit, however, 
cannot be taken from the finite-e expressions we have derived for particular profiles, since their derivation 
assumes that, for all x^., cr^^Xd) > 0, in order to map R to a finite interval, finding the eigenfunctions in that 
finite region. The use of that mapping to find those eigenfunctions collapses when e — > 0, however, since 
then aJxd) becomes a (5- function. 



B. Finite width mirrors 

We now study some general properties of £o{a), which are independent of any assumption regarding the 
form of the defect. 

We first note that, denoting by Ow = {^l!l\0\'4'l') the matrix elements of the operator O in the basis of 
eigenfunctions of IC, and a similar convention for the V and Q kernels, we have: 

Ow = — VP/™ — Qrni' . (97) 

ai ^ am 

To study the properties of Sq in more detail, we consider the explicit form of Vim and Qim- 
Vim ^Jd^J dz'^i{z)e-^\^-(^-^-^-(^'^-'^\^Pm{z') 

Qim = ^jdzj dz>K2)e-'^'"''("'-"''^^'^+"l^m(^') = Vml • (98) 

0„, = 1 y . (99) 

ai ^ am 

ra 

An eigenvalue -d corresponding to O (required in order to evaluate the Casimir energy by applying (j50[) ) is 
then determined by Owvii ~ -dvi. Using (|99p . and taking into account the fact that ai > 0, this may be 
equivalently written in matrix form as follows: 

Av = dgv (100) 



Then, 
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where A is a symmetric and positive definite matrix whose elements are given by 

Am = (101) 

and g is 'A diagonal matrix: g = diagjai, Q!2, • ■ •}• 

Note that (|100p has the form of a generalized eigenvalue problem for the symmetric matrix A and the 
diagonal positive matrix g. Taking advantage of the fact that g is diagonal, one can show that the eigenvalues 

of A can be found as the ones of a matrix B such that Bki = af.'^ Akia^ ^ . This matrix B may be given an 
even more convenient form: 

B = CC\ (102) 

where 

Cm = . (103) 

Being symmetric and definite positive, there are many available numerical algorithms to compute the eigen- 
values of the matrix B. 



V. CONCLUSIONS 



In this paper we have computed the Casimir energy for a real scalar field in different backgrounds that 
describe finite width, semitransparent mirrors. The properties of the mirrors are described by A(a;o — 
Xq;X|| — Xy) and CTe(xd). The former is a two-point function that depends on time and on the parallel 
coordinates, and represents frequency-dependent transmission and reflection coefficients of the mirrors, while 
the latter depends only on the coordinate normal to the mirrors, and describes the spatial dependence of 
their electromagnetic properties. Our starting point, a non local effective action for the quantum scalar field, 
can be thought as arising from the interaction with the degrees of freedom of the mirrors. 

The functions A and fTg act as physical regulators for the divergences of the zero point energy, and the 
usual 'perfect conductor' limit can be obtained when becomes a (S-function and A tends to infinity. In this 
case, the effect of the mirrors is to impose Dirichlet boundary conditions on the scalar field. 

We have computed explicitly the self-energy of a mirror using smooth and piecewise constant profiles a^^. 
In both cases, the self-energy is finite under the assumption that the mirror becomes transparent at high 
frequencies. We have discussed in detail the dependence of the results with the cutoff frequency Wc at which 
the mirror becomes transparent, and with the width e of the mirror, showing that the zero point energy 
vanishes for eujc < 1. Finally, we have also analyzed the UV divergences that arise as e ^ 0, recovering 
previous results Q . We hope that an analogous computation of local quantities like < T^^, > will, in this 
context, be useful to discuss the gravitational effects of the zero point fluctuations, without having to deal 
neither with bulk nor surface divergences. 

For the case of two mirrors, we have discussed some general properties of the interaction energy, for 
arbitrary functions A and a^. In the thin wall limit, the expression for the energy is considerably simpler, 
and we have evaluated it explicitly for some particular cases, reproducing also the well known results for 
perfect mirrors. 

Regarding future research, as already mentioned, it would certainly be of interest to compute also local 
quantities, like the energy and pressure densities. Moreover, a realistic model for the interaction between the 
quantum field and the degrees of freedom in the mirror would allow us to derive a nonlocal effective action 
suitable for a detailed analysis of the dissipation effects. This realistic model would necessarily involve the 
electromagnetic field. Due to gauge invariance, on general grounds we expect the nonlocal effective action 
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to be of the form: 

5eff = j d'^+^x j d^'+^x' F^,(x)if^^'"^(x;x')i^pa(x') , (104) 

where the kernel K'^'^'"^ {x; x') encodes the electromagnetic properties of the mirrors. Here gauge invariance 
is inherited from the (assumed) gauge invariant coupling between the material and the gauge field. On the 
other hand, this is consistent with the limiting case of ideal plates, where the boundary conditions are given 
in terms of F^i, rather than . In spite of the fact that the effective action shall involve derivatives of the 
gauge field, we expect a conveniently adapted version of the method will make it possible to consider also 
this case. 

The extension of the results of this paper to non planar and/or non static mirrors would also be of high 
interest. 
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